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to the streaming direction is the fastest growing mode, and probably the first to appear, consistent 
with the aligned filaments that are seen in PIC simulations. Oblique, quasiperpendicular modes 
grow almost as fast, as the growth rate varies only moderately with angle, and they may distort or 
corrugate the filaments after the perpendicular mode saturates. 
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I. INTRODUCTION 



The Weibel instability appears to be important for catalyzing shocks in unmagnetized plasmas. With the discovery 
of relativistic blast waves from GRB, it has received much attention, [l|-[8| , and appears to be confirmed by simulations. 
The conventional wisdom is that when two oppositely directed plasma streams collide, the Weibel instability generates 
small-scale magnetic fields; particle scattering off these magnetic fluctuations provides an isotropization mechanism 
necessary for the shock transition to form. Particle in cell simulations |9l.ll0j show that a) magnetic filaments grow with 
a k vector that is perpendicular to the beam direction, b) these filaments saturate at about 10% of the equipartition 
level and are eventually disrupted, and c) the electrons are efficiently heated within the shock structure until their 
energy becomes comparable to the proton energy. The full shock transition occurs at a small scale of dozens proton 
skin depths. The physical reason for the remarkably efficient electron heating has not been fully established yet, 
however the transfer of energy to the electrons can already be seen in the linear instabilities, which are accompanied 
by induced electric fields that induce electron countercurrents, as per Lenz's law 0,0, EH- ^ is necessary to understand 
the underlying physics before simulations can be scaled and generalized to real astrophysical phenomena. 

An important point regarding the very first stages of relaxation is that the electron streaming is halted easily, 
whereas protons still plow on through an isotropic electron gas. The newly isotropized electron background effectively 
suppresses the generation of magnetic field because of the induced electric field and attendant Lenz currents opposite 
to the proton currents [!, 0, ml In the case of a Weibel unstable mode with k vector exactly perpendicular to 
counterstreaming ion beams, the mode has a purely imaginary frequency in the electron frame. The magnitude of 
the Lenz currents then decreases with electron temperature. In non-relativistic case, the fraction of electrons that 
resonate with the purely growing mode decreases with the increasing thermal velocity. In relativistically hot plasma, 
the electron inertia grows with the plasma temperature. In any case, both the magnitude and the spatial scale of 
the magnetic fluctuations depend on the "temperature" of the background electrons; the larger the temperature, the 
stronger the magnetic fluctuations. 

As a step towards understanding the filament disruption and electron heating within the shock structure, one 
should consider the full spectrum of unstable modes in the presence of proton beams within the relativistically hot 
electron background and calculate how much of the energy is transferred from protons to electrons. Two electron 
heating mechanisms are possible in the linear stages of instability: 1) Lenz currents, mentioned above, and 2) Landau 
damping against hot electrons that move in resonance with the wave, feeling its longitudinal electric field along the 
direction of the wave's phase velocity vector. The second happens only when the electrons are already hot enough to 
be in resonance with a wave of substantial phase velocity. These are the waves that make an oblique angle with the 
beam direction. The first can occur even for cold electrons, because the phase velocity of the perpendicular Weibel 
unstable mode can vanish in the frame of the electrons. Thus the relative importance of the two mechanisms hinges 
on the role of oblique waves, and the effects of hot electrons on them. 

It has been claimed recently that the instability of oblique modes [H and the longitudinal Buneman instability 
compete with the Weibel instability in the relativistic case. In this paper we study the full spectrum of the unstable 
modes for a relativistic beam in relativistically hot (T ^> m e c 2 ) electron background, and we show that the electrons 
Landau damp the oblique modes. The unstable oblique modes then take on more of Weibcl-like character, with 
corresponding growth rate. 

The case of oblique modes has already been studied long ago for cold plasmas. Here, specific properties of relativistic 
streaming instabilities are dictated by the fact that the response of the relativistic particles to an external perturbation 
is highly anisotropic. The effective particle mass, namely the ratio of force to acceleration, for longitudinal (with 
respect to the direction of the particle motion) perturbations is m c ff = m^ 3 whereas for the transverse perturbations 
m e ff = 7717; here 7 is the particle Lorentz factor. For this reason, a relativistic particle beam excites, in a cold plasma, 
obliquely propagating Langmuir waves much faster than the waves in the direction of the beam (l2l . [l3| : the growth 
rate of the resonant beam instability is F oc m cff 1 ^ 3 so that F cx 7 -1 for longitudinal modes and T oc 7~ 1 /3 for oblique 
modes, respectively. Fully electromagnetic analysis of a cold plasma and a cold beam system FJ| confirmed that the 
obliquely propagating unstable modes are practically electrostatic and also revealed an instability of purely transverse 
perturbations with the growth rate T oc m cff 1 ^ 2 oc 7 -1 / 2 , which is in fact the relativistic counterpart of the Weibel 
instability (l5j . 

The fact that a relativistic beam is most unstable to the excitation of oblique waves has important implications 
for inertial fusion (e.g. [l6[), where the electrons remain subrelativistic. In a relativistically hot plasma, however, 
electrons can Landau damp oblique waves, which travel at a significant fraction of c, and the question should be 
investigated carefully. 

Another interesting aspect of the simulation results is that current filaments aligned with the beams are formed 
initially, but eventually disrupt. The disruption may be due to non-linear effects, but the filaments may also be 
disrupted by the eventually emergence of oblique modes, which one suspects would lead to a more complicated 
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magnetic geometry and chaotic particle trajectories. The randomization of electron trajectories is of course connected 
to the issue of electron heating, as the rate of heating by electric fields is determined by electric resistivity, which is 
proportional to the electron scattering rate. So oblique modes may play a role both in electron heating and filament 
disruption. 

In ref. Q, the general dispersion equation was solved numerically and the results were reported graphically for 
a specific sample choices of parameters. Here we use a more analytic approach, in the interest of generality. We 
present analytic solutions to the general dispersion equation describing a monochromatic beam in the relativistically 
hot electron background. Generalization to multi-beam systems is straightforward. The mass of the beam particles 
is not specified here, but is intended in the context of proton beams within the shock structure. 

It is shown that the effects of Landau damping are comparable to or less than those of induced transverse electric 
fields, so that the general picture that is gaining acceptance - beam instability that is primarily Weibel in charac- 
ter, generating magnetic filaments of negligible phase velocity - more or less holds up even when oblique modes and 
Landau damping are included in the analysis. On the other hand, the oblique modes, as they become quasiperpen- 
dicular, increasingly resemble Weibel unstable electromagnetic modes, and have a growth rate that is a bit slower but 
comparable to the latter. So they are not entirely negligible and could eventually lead to a more complex magnetic 
structure. 

We consider here hot electrons but continue to assume cold ion beams. When the ion beams are warm, all instabilities 
are suppressed, but the Weibel branch is still unstable at sufficiently long wavelengths, where the electrostatic modes 
are superluminal and stable. So we consider the cold ion beam case the most favorable for the electrostatic instabilities. 
We will see that even in this case, the Weibel instability is dominated in the relativistically hot plasma. 

The paper is organized as follows. In the next section, we present the dispersion equations for the monochrmatic 
beam in the isotropic electron background. In sect. Ill we study instability of longitudinal (along the beam direction) 
modes, the Buneman instability. The purely transverse mode corresponding to the Weibel instability are considered 
in sect. IV. Stability of oblique modes is addressed in sect. V. Conclusions are presented in sect. VI. In Appendix, 
we present the permeabilities of the electron plasma with ultrarelativistic Maxwell's distribution. 



II. DISPERSION EQUATION 

We consider interaction between a relativistic "monochromatic" beam and the isotropic electron plasma. The 
stability analysis is reduced to solution of the dispersion equation 



Det|/c 2 <V 



kikj to 



o. 



(i) 



In this paper, we take the speed of light to be unity. In the plasma-beam system, one can conveniently separate the 
dielectric tensor into the plasma and the beam parts, Sij = + £y — <V,-, and present the plasma dielectric tensor as 
(e.g. 0) 



4 = ( 5 »- 



k<i k j 



(2) 



where et and ei are the transverse and longitudinal dielectric permeabilities, respectively. The permeabilities of the 
Maxwell plasma with highly relativistic temperatures are presented in Appendix. We take a monochromatic beam of 
particles with the mass rrif, and the density rib moving in z direction with the velocity Vb- Then dielectric tensor of 
the beam is presented as 



4 = M 1 



w 2 7b 



^b ( k i V bj + fyVbi 



U 2 Jb 



u - k z v b [u ~ k z v) 



-VbiVbj 



(3) 



where 



uj\ = iirnbe 2 /nib- 



(4) 



We choose the coordinate system such that k y = 0. Then the dispersion equation for the plasma-beam system is 
reduced, after some algebra, to the form 



7^ (cj - k z v b y 



u 2 e t + -± 



7 (lu - k z v b y 



(5) 
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If the plasma is cold, e t — £; = 1 — uj 2 /uj 2 , uj 2 = Aire 2 n/m e , equation ([5]) reproduces the result of 
For all practical purposes we can safely assume that 

" b «-^=. (6) 



For electron beams, this implies Tif,/% <C n p /^T whereas for proton beams, which is of special interest in the context 
of relativistic shocks, this condition is fulfilled even for n b ~ n p , 7 & ~ j T . 

III. BUNEMAN INSTABILITY 

First we consider instability of the wave propagating along the beam, k = (0,0, k). In this case the dispersion 
equation §5§ is reduced to 

ei = y. (7) 

% (w - 

The instability occurs due to the resonance of particles with waves, so that one can write w = kvb + 8u>; 8lj <C w. 
For a highly relativistic beam, w ss fc, therefore we have to use equation pip for e;. Then the dispersion equation is 
reduced to the cubic equation 



Slu 



Slu + £(fc - fc ) - 7^3 



fco 

H 



^b^O /OS 



The maximal growth rate is achieved at fc = ko[l + l/(2^ 2 )] and is equal to 

3V«/ W6 \' Wo 3 1 / 6 ^ 2/3 ^y 3 VhT27¥ 
BU ~ " 2V3 U J 7^ /3_ 25/3 7,7^ ' 

The Buneman regime occurs at the condition ^ > 7^. When the beam velocity is less than the plasma thermal 
velocities, jb *C 75™, the beam is still unstable, but the growth rate decreases. According to ref. [19j . the maximal 
growth rate in this case is estimated as 



r « (10) 

7 fc 3/2 



IV. WEIBEL INSTABILITY 



Now let us consider instability of transverse waves, k = (fc, 0,0). In this case one can neglect the terms with LUb 
in the left-hand side of the dispersion equation (JS|). The right-hand side becomes comparable with the left-hand side 
only at w <C k. In this limit, the plasma dielectric permeabilities (|27p and ([25)1 are reduced to 

2 99 
1 + ■ 1 + ' ^ (11) 

Now the dispersion equation (|S|) is reduced to the cubic equation 

1 + — ^ fc 2 - = 12 

V 7tk y V 47t/c / 7T7b^ 

The unstable solutions to this equation are purely imaginary, which implies aperiodic instability. The growth rate 
has a broad maximum (see FigfT]) 

Fweibel ~ —= (13) 
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FIG. 1. Growth rate of the Weibel instability of the proton beam at 75 = 100, 7t = 20 and equal densities rib = n 



at 



2 \ 1 / 3 



1/2 



(14) 



In the long and small wavelength limits the growth rate is given by: 



,2„,_\ 1/3 




Comparing the Buneman and the Weibel instabilities, one finds that the ratio of the growth rates very weakly 
depends on the mass and density ratios of the beam and plasma particles: 

3 1 / 6 foj p \ 1/3 VhT277 _ 3^6 fm b n p \ 1/6 0n2^ 



Tweibel 25/3 ^ J ^5/6 \m e n b J J^ 2 ^ 6 ' 

Therefore in the relativistic case, 7t,7& 3> 1, the Weibel instability dominates the Buneman one even for the proton 
beam unless the density of the beam is extraordinarily small. 



V. INSTABILITY OF OBLIQUE MODES 

Now let us consider oblique modes. In this case one can also neglect the terms with lui, in the left-hand side of the 
dispersion equation ([5])- The right-hand side becomes comparable with the left-hand side only for modes in resonance 
with the beam, i.e. if u) ~ k z = kcosO. Therefore one can substitute uj = kcosO into all the terms with the exception 
of the resonance denominator, which yields 

uj 2 sin 2 

ei(i - cos 2 ee t ) = - 2 1 1 _ sin2 9e i - co&2 M > ( 18 ) 

7 (w — kvi, cos 6) 



FIG. 2. Growth rate of the oblique modes as a function of k x ,k z . The beam is in the z direction. 



where one can use the permeabilities in the form (|2"T)) and (|28p with w = k cos 9: 



UJ 2 ( 9 7T \ 

Ei = l + -r7r— l-cos01ncot- + -icos0 (19) 
k z jr V 2 2/ 

et = 1 f sin 2 6» In cot - + cos - -i sin 2 0^ . (20) 

2fc 2 7T cos0 V 2 2 / v ; 

Note that these expression diverge at 9 = and 6* = 7r/2; therefore these specific cases should be considered separately; 
this has been done in sect. Ill and IV. 

Now the solution to the dispersion equation is written as 

u = kcos9± -^=$ [ -^,0 J ; (21) 



76 \ irk 2 



' 1 — sin 2 6 Ei — cos 2 9e t 



■ ... — COS Utt , n 

$ = sm^ n 2/3 \ ■ 22 

1 £/(l — cos 9e t ) 



The growth rate of the instability is presented as 



2 



V76 \ 7Tfe 2 / 



One sees that the dependence on the beam parameters is split from the dependence on the plasma parameters therefore 
we can find the universal dependence of the growth rate on the parameters by normalizing T by Uf,/ and k by 
LU p/\/lT- This dependence is presented in Fig. ([2]). The polarization of the unstable modes is hsown in Fig. ([3]). In 
Fig. ((U), the growth rate is presented as a function of k at a fixed angle 9 = 7r/4. In Fig. ([S]), we plot the maximal 
growth rate as a function of the angle between the wave and the beam. One sees that the largest growth rate is 
achieved for the transverse mode, i.e. for the Weibel instability. 



VI. DISCUSSION AND CONCLUSIONS 



We found the full set of unstable modes for the monochromatic relativistic beam propagating through relativistically 
hot electron background. The results are easily generalized to the multi-stream case. The oblique and longitudinal 
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FIG. 3. The direction of the electric field in the oblique modes; 76 = f00, 7t = 20, n = 2nb- Note that as k decreases, the 
modes become more electromagnetic in character. 
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FIG. 4. Growth rate of the oblique mode as a function of k at 8 — 7r/4 



modes are resonant, therefore each beam excites the appropriate wave independently of others. In the case of the 
Weibel instability, one has to substitute, for the term uif/^b in the right-hand side of the dispersion equation (fl"2j).the 
sum over all the beams ■ u>? Jjb j] this yields the same substitution in the expressions for the growth rate (see also 

Our analysis confirms that the strongest is the Weibel instability with wavevector perpendicular to the beam 
direction, with the growth rate given by equation (|13|) . The reason the oblique modes are suppressed, relative to the 
case of a relativistic beam in a non-relativistic plasma, is that the resonance waves are subluminal, oj/k = vi,cos9, 
therefore in a relativistically hot plasma, they are suppressed by the Landau damping. The reason, on the other 
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FIG. 5. The maximal growth rate as a function of 9 



hand, the oblique modes grow more slowly than the perpendicular Weibel instability, which is also suppressed by the 
electrons, is that the projection of the beam anisotropy on the k vector is reduced by sin0, so such suppression is 
modest when 9 is close to ir/2. 

Support from the Israel Science Foundation, the Israel-U.S. Binational Science Foundation, and the Joan and Robert 
Arnow Chair of Theoretical Astrophysics is gratefully acknowledged. 



APPENDIX. DIELECTRIC PERMEABILITIES OF RELATIVISTICALLY HOT PLASMA. 



Generally the longitudinal and transverse permeabilities are written as (e.g. [13]) 



sf = l 



-pi _ 



1 + 



Aire 2 



Aire 1 
k- v 



df(p) 



d 3 P 



k • 



dp uj + iO — k • v 

0/(p) d 3 p 
dp uj + iO — k • v ' 



We assume that the plasma electrons have Maxwell's distribution 

71 

/(P) 



87r(m e 7T) 3 



exp 



J_ 

It 



(24) 
(25) 

(26) 



with ultrarclativistic temperatures, 3> 1- Integration over the angles in Eqs. (|24|) and ([25]) is performed straight- 
forwardly. Integration over p could be performed if one substitutes v = 1; this yields [18( (see also [171 ]) 



1 + — In 

2k 



uj + iO — k 
uj + k 



sf = l 



iujkjq 



uj + iO — k 
uj + k 



2uj 
~~k 



where 



4ire 2 n 



uj; 



(27) 
(28) 

(29) 
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Note that the expression (|27| diverges at u> —> k. In this case, one has to take into account small deviations of v 
from the speed of light in denominators u> — k • v in the exact formula (|24p . Only then one can find correctly describe 
the longitudinal waves with the phase velocity close to the speed of light and in particular, subluminal modes crucially 
important for the beam instability [l9j . Following [l9[ we first find the longitudinal wave with the phase velocity 
equal to speed of light uj /k — 1. In this case the dispersion equation, si — 0, with the exact permeability (f24f yields 



w5 = fcg = — ln2 7T . 

It 



(30) 



Now we could expand si in small w — uio, k — ko to yield the longitudinal permeability in the vicinity of the point 
cj ,fc : 



ef(u, k) 



■ 1 (w - Ul ) + 1 ' 



dui 



dk 



(k - k ). 



(31) 



o \ / o 

The zero subscript means that the derivatives are taken at k = feo, uj = loq. Straightforward calculation yields 



'def\ _ 1%4jt 

6U> I Wn 



(32) 



and 



can be represented as 



dk 



= -(i-0 



dui 



(33) 



where £ is some function of temperature, it is small as 7 T 2 (see Fig. [5]). Since ( ) > 



dk 



the waves with 



k > ko are subluminal; these waves could be resonantly excited by a relativistic particle beam, see sect. III. 
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FIG. 6. The parameter £ (see Eq. (|33[) 1 as a function of the plasma temperature 



